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I. INTRODUCTION. 

Linear in temperature corrections in the spin suscepti- 
bility of the two-dimensional (2D) electron gas has been 
discussed intensively in the past decade.— &MiMiZJL1 Our 
interest to this question is motivated by a recent observa- 
tion of a strong temperature dependence in the spin sus- 
ceptibility in the silicon mctal-oxidc-scmiconductor ficld- 
effect transistor (Si-MOSFET)3 It has been noticed by 
Misawa^ that linear in T corrections, if exist, should be a 
result of a non-analytic behavior of the thermodynamic 
potential at a small temperature and magnetic field. In- 
deed, a naive Taylor expansion term H 2 T in the thermo- 
dynamic potential corresponding to the linear in T cor- 
rection to the spin susceptibility would violate the third 
Law of thermodynamics. Therefore, there should be a 
strong dependence on the order of limits H — ► 0, T — > 
in the thermodynamic potential. Obviously, to analyze 
such a non-analytic behavior one has to go beyond the 
standard theory of the Fermi liquid systems. Still, we 
use a machinery of the microscopic Fermi liquid theory 
as our starting point. 

It is known that a pair of Green's functions has a sin- 
gular behavior when it has a momentum transfer close 
to 2pp. In the polarization operator this singularity re- 
veals itself as the Kohn anomaly at 2pp. In the case of 
2D, the Kohn anomaly leads to a cusp- like dependence 
on the transferred momentum and frequenc y 11 ' 12 ! 13 This 
cusp-like behavior is related to the fact that for a 2pp 
momentum transfer the two patches on the Fermi sur- 
face are located on the opposite sides and are parallel 
to each other. The cusp-edges are sensitive to a relative 
shift of the chemical potentials. When the Fermi sur- 
face is spin-split due to magnetic field, the sharp effect 
of shifting of the two cusps near the point when they 
touch generates a non-analytic temperature dependen- 
cies in the thermodynamics quantities. This touching of 
singularities gives a clue for understanding why there is 
an anomaly in the spin-susceptibility. The analysis of the 
temperature correction with the use of the non-analytic 
parts in the product of two pairs of Green's functions 



discussed in this paper is different from the previous cal- 
culations; e.g. Ref. @. We believe that it makes the 
origin of the non-analytical behavior of the thermody- 
namic potential more transparent. The dependence on 
the order of limits H, T — > are explicitly demonstrated 
in this calculation. 

The calculation of the non- analytic term in the thermo- 
dynamic potential containing two pairs of Green's func- 
tions shows that the anomalous temperature dependence 
originates from small frequencies and momenta close to 
the Fermi surface. The momenta are restricted to the 
narrow angular intervals corresponding to scattering with 
a 2pp momentum transfer as it is shown in Fig. [TJ The 
amplitude associated with such scattering processes is 
natural to be called the backward scattering amplitude 
and will be denoted as r(7r). Because of a special config- 
uration of its four momenta, the amplitude r(7r) can act 
also within the Cooper channel^ Moreover, the product 
of two 2pp-pairs of Green's functions by an appropriate 
twisting of the Greens functions can be read as the prod- 
uct of two pairs in the Cooper channel (or two particle- 
hole pairs in the zero-sound channel^). In this paper 
we demonstrate that the non-analytic terms, associated 
previously only with 2pj?-scattering, should be analyzed 
having in mind the rescattering of pairs in the Cooper 
channel^ We emphasize here the rescattering in the 
Cooper channel as the only source of the temperature- 
dependent renormalizations of the thermodynamic quan- 
tities of the Fermi liquid systems in the ballistic regime. 
The analysis of the renormalization of the anomalous 
temperature dependencies is performed in terms of an- 
gular harmonics in the Cooper channel, because each 
harmonic is renormalized differently. In the calculation 
of the linear in T term with the use of harmonics in the 
Cooper channel we have not assumed a priory the special 
importance of the backward scattering amplitude r(?r). 
Interestingly, we obtain this fact as a result of indepen- 
dent calculation of the product of two pairs. 

We show that in the Cooper channel ladder the non- 
analytic terms in the spin-susceptibility, 5%, are gener- 
ated by the products with arbitrary number of the pairs 
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of Green's functions, i.e., not only a product of the two 
pairs. This happens because a correlation function de- 
scribing the propagation of a pair of quasiparticles in the 
Cooper channel has a dynamic part which depends on the 
ratio uj/qvF- (In the microscopic Fermi liquid theory^ 
the scattering amplitude in the zero sound channel has 
exactly the same feature.) Consequently, we obtain that 
dSx/dT is given by a power series in the renormalizcd 
Cooper channel amplitudes that does not reduce to the 
renormalization group generalization of the result ob- 
tained in the second order. The truncation of this series 
is possible only when the Cooper channel amplitudes are 
small. 

This paper is organized as follows. In Section 2 a 
general discussion of the anomalous temperature depen- 
dences in the Fermi liquid systems is presented. We start 
by stressing that to get anomalous temperature depen- 
dencies it is necessary to pinpoint non-analyticities in 
the thermodynamic potential to prevent regular Som- 
merfeld's expansion at low temperatures. To get a clue 
of the origin of the non-analytic terms, an analogy with 
an auxiliary model system is presented. In the end of 
this Section wc explain why in the analysis of the non- 
analytic terms in the thermodynamic potential the lad- 
der diagrams are of particular importance. In Section 3 
the physical consequences of the renormalizations in the 
Cooper channel are discussed in connection with the tem- 
perature dependence in the spin susceptibility observed 
in the Si-MOSFET^ It is argued that at low tempera- 
ture when the interaction amplitudes arc strongly renor- 
malized their own temperature dependence may over- 
come the linear in temperature factor. This may be a 
possible explanation of the observed sign of the tempera- 
ture dependence of the spin susceptibility in the 2D elec- 
tron gas. All technical details of the calculations are 
moved in Appendices. We first present details of the 
calculations of the linear in T term in the spin suscepti- 
bility originating from two 2pi?-pairs of quasiparticles. In 
Appendix [B] we obtain an expression which demonstrate 
explicitly how the dependence on the order of taking the 
limits H — ► 0, T — > appears in the thermodynamic 
potential (for more discussions on this subject see also 
Ref. [13 )■ Next, in Appendices [Cl and iDl the tempera- 
ture corrections arc analyzed within the Cooper channel. 
In the concluding Section 4 we give several historical re- 
marks. In particular we stress a rather unique nature of 
the temperature dependences in 2D. 



II. GENERAL DISCUSSION OF ANOMALOUS 
TERMS IN x- 

Let us discuss the origin of the anomalous temperature 
dependences. The calculation of interaction corrections 
to a thermodynamic quantity requires a summation over 
bosonic frequency iu> n . Transforming the sum into an in- 
tegral along the real frequency axis results in an integral 
of the form J dco coth(/?w/2)/(a;). If f(u>) is smooth and 




FIG. 1: The diagrams with four Green's functions in the ther- 
modynamic potential which produce the anomalous term in 
the spin susceptibility. p+ and p_ indicate the directions of 
the momenta of electrons: p+ « — p_. The chemical poten- 
tial shifts and S/j.2 has been introduced in each of the two 
pairs of Green's function. The spin-susceptibility can be gen- 
erated from these diagrams by taking the second derivative 
with respect to <5^i,2- The triangles attached to the external 
vertices indicate the Fermi-liquid renormalization of the static 
spin vertices. 



regular in the vicinity of w = 0, the standard Sommerfcld 
expansion will involve only even powers of temperature. 
In order to get a term with an odd power of temperature 
(e.g., linear in T correction to the spin-susceptibility), / 
should be non-analytic preventing the Taylor expansion 
at lo = 0. Let us see how such non-analyticity develops 
in the case of the 2D electron gas. It is shown in Appen- 
dices [B] and [D] that both a 2pp--pair of Green's functions 
and a pair in the Cooper channel have non-analytic parts 
(edge parts). Together these edge parts act as a sort of 
an "edge-mode" propagating in the opposite directions. 
In the case of a 2pj?-pair they are denoted as II r and 
II;, i.e., "right" and "left", reflecting the combinations 
lo ± VF^q in the denominators of Eq. (|A9[) . In the case 
of a pair in the Cooper channel the left and right edges 
of the branch-cuts in if? play a similar role, see Eqs. (|D2|) 
and (|D10p . A diagram in the thermodynamic potential 
with the two edges running on each other produces the 
anomalous temperature dependence in the spin suscepti- 
bility because of its high sensitivity to the temperature 
and magnetic field at the point when they touch, see 
Fig. [21 The spectral weight of the edge-modes is much 
weaker compared to that of the usual quasiparticles, e.g., 
phonons. However, taking derivative with respect to an 
external parameter, a magnetic field in the case of the 
spin susceptibility, makes the spectral weight singular. 
As a result, the edge modes become active and are able 
to generate a non- analytic / at lo = 0. 

To get a clue of the origin of the discussed non- 
analytic terms, the following analogy may be useful. As 
it has been mentioned, a product of two propagators of 
the edge modes running on each other generates non- 
analytic terms in the spin susceptibility. With this in 
mind, consider an ensemble of one-dimensional right- and 
left-moving particles turning into each other due to a 
mixing matrix element. In such an auxiliary system, 
the diagrammatic expansion with respect to the mix- 
ing amplitude would suffer from strong "infrared" di- 
vergences due to the degeneracy of the right- and left- 
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moving particles at q = 0. In particular, in the second 
order in mixing the thermodynamic potential has a term 
~ Vmix T T,J M^l + 9 2 )" 1 - The avoiding of the level 
crossing leads to the restructuring of the ground state of 
this system with the gap opening at q = 0. Coming back 
to the non-analytic term in spin susceptibility in 2D at a 
finite temperature, we can see that the magnetic field H 
(multiplied by the interaction amplitude) acts as the mix- 
ing matrix element in the subsystem of the edge modes. 
The considered diagrams are not divergent, because the 
spectral weight of the edge-modes is much weaker com- 
pared to the usual quasiparticles, e.g., phonons. There- 
fore, there are no enough resources for restructuring of 
the ground state. In addition, being composite objects, 
the edge modes are smeared by the temperature through 
the Fermi-Dirac function np of the two fermions from 
which they are made of. Nevertheless, they still are able 
to produce the anomalous terms in the thermodynamic 
potential which appear as a trace of the mentioned above 
divergences in the expansion with respect to mixing am- 
plitude in the auxiliary system. The analogy with auxil- 
iary system also helps to understand the dependence on 
the order of taking the limits H — > 0, T — > 0. This de- 
pendence is similar to the competition between the value 
of the mixing matrix element and the energy of the par- 
ticles in the auxiliary system, provided that the energy 
is substituted by the temperature in the original system 
of interacting electrons. 

In Apppcndices we show how the non-analytic (edge) 
parts in the product of two pairs of Green's functions 
generate anomalous dependencies in the thermodynamic 
potential. With this experience let us now discuss possi- 
ble generalizations of this mechanism. In diagrams defin- 
ing the spin susceptibility there are two differentiations 
with respect to the magnetic field. The calculation of 
the term with two pairs indicates that for nucleating 
the non-analyticity in the thermodynamic potential the 
quasiparticles involved in the differentiation have to be 
constrained by the conservation of momentum and en- 
ergy. For the product of two pairs these constraints are 
imposed automatically. Since for the more complicated 
diagrams this may not be the case, in identifying other 
non-analytic terms it is important to preserve this fea- 
ture. Consequently, in the thermodynamic potential only 
those diagrams are essential that can be arranged as a 
closed ladder loop, and where only the Green's function 
within the sections of the ladder are differentiated with 
respect to the magnetic field (by section we understand 
a pair of Green's function describing a propagation of a 
pair of quasiparticles between the scattering events). 

Organizing a diagram in the form of a ladder can be 
performed as follows. (In the discussion below the closed 
ladder loop should be not confused with the fermionic 
loops). If the Green functions marked by the differenti- 
ation belong to two different fermionic loops and there 
exists at least one another pair such that removing this 
group of Green's functions (articulation quadruplet in the 
terminology of Ref.fisl) will split the diagram into two dis- 



connected pieces, then these Green's functions nucleate 
the ladder. Obviously, the Green's functions that split 
the diagram are constrained by the conservation of en- 
ergy and momentum. When such splitting requires more 
than one supporting pair, the non-analytic parts, which 
are a goal of our studies, are smeared by the additional 
integrations and this diagram cannot generate the lin- 
ear in temperature term in spin susceptibility. When the 
two differentiated Green's functions belong to the same 
fermionic loop, one first chooses another fermionic loop 
and repeats the procedure. Namely, it is needed to find 
another supporting pair of Green's functions within the 
second loop which allows splitting the diagram into dis- 
connected pieces. This completes the procedure of ar- 
ranging the differentiated Green's functions as a part of 
a ladder. 

Unlike the one-dimensional case, the renormalizations 
in the Cooper and 2pp channels are not mutually cross- 
coupled by the logarithmic corrections in two dimensions. 
This makes the Cooper channel the only one in which 
terms logarithmically divergent in temperature are gen- 
erated. Besides that, in the treatment of the Cooper 
ladder we can consider its crossbars while ignoring their 
non-analytic parts originating from the 2pp channel. We 
rely here on that (i) the non-analytic parts in the 2pp- 
scattering are small unless they are differentiated (notice 
that all the differentiated Green's functions have been 
attributed to the Cooper ladder here), and more impor- 
tantly that (ii) the non-analyticity generated in another 
channel is smeared out. That's why for the calculation 
of the anomalous terms generated in the Cooper chan- 
nel it is enough to consider a ladder assuming that the 
crossbars in the Cooper ladder are analytic— On these 
grounds, for the analysis of the non-analytic temperature 
dependences we can use as a starting point the conven- 
tional theory of the Fermi-liquid systems. 



III. PHYSICAL CONSEQUENCES OF THE 
RENORMALIZATIONS IN THE COOPER 
CHANNEL. 

Cooper's logarithms lead to renormalizations of the lin- 
ear in T term in the spin-susceptibility generated in the 
Cooper and 2pp channels (see Appendix [Cl for details): 

<5x = ^-|r c (7r) T | 2 , (1) 

where v = m/2ir is the density of states per one species 
(to is the effective mass which includes Fermi liquid 
renormalizations), and we omit the factors (gp,B/2) ev- 
erywhere. The renormalized backward scattering ampli- 
tude r c (7r) is 

^)r = D-M. ^= l + rcL/T - 

(2) 
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Here "f n (T) are the renormalizcd Fourier harmonics am- 
plitudes in the Cooper channel; r„ are bare amplitudes 
defined at a large energy scale ~ cf ■ This result is in full 
correspondence with the renormalization group equation 
for the angle-dependent amplitude T c (0) in the Cooper- 
channel. The fact that the amplitudes 7n(T) undergo 
strong rcnormalizations with temperature leads to im- 
portant physical consequences which we now discuss. 

When the temperature is reduced, the amplitudes 
7 n (T) which have a positive value are renormalized to 
zero, while those that initially are negative grow in mag- 
nitude. The initial (bare) value of T c ' (-k) is most prob- 
ably dominated by the zero-harmonic amplitude. It is 
naturally to expect that this harmonic is repulsive, and 
therefore when the temperature is lowered it dies out. In 
the intermediate region of temperatures when negative 
amplitudes are still small, i.e., when l/|r„| > InA/T > 
1/r^, the temperature dependence of the spin suscepti- 
bility is 

, , . T 



This result has been also obtained in Ref. [2(| using a 
technique developed recently for the ballistic systems— 
As the temperature is lowered, the negative amplitudes 
start to grow. There are general reasons why the negative 
amplitudes are always present for some harmonics— The 
largest negative amplitude, j n = 7 C , is most important 
as this amplitude grows most rapidly approaching the 
Cooper instability. When the amplitude |7 C (T)| becomes 
of order of unity, it acquires an essential temperature 
dependence. At small enough temperatures, when the 
renormalized to zero positive harmonics are suppressed, 
see Eq. Q, the temperature dependence of the grow- 
ing negative amplitude j c (T) = T c /(1 + T c InA/T) 
prevails and can determine the temperature behavior of 
the spin susceptibility. This may cause a non-monotonic 
behavior in 5\{T). For the purpose of illustration, we 
will discuss the temperature dependence in the spin sus- 
ceptibility leaving only the most rapidly growing ampli- 
tude 7 (T). In this case the temperature dependence 
of Sx(T) oc Tj c (T) 2 is determined by two competing 
terms: 

dS X /dT = (u/e F )[ 1 c (T) 2 + 2T 1 c (T)(d 1 c (T)/dT)} 
= ( V /e F )( 1 c (Tf~2\ 1 c (T)f). (4) 

This expression changes its sign when — j c (T) > 1/2. 

Let us now consider the present experimental situa- 
tion^ The data indicate a noticeable temperature de- 
pendence of the spin susceptibility at T > 2K, which 
is too strong to be attributed to the conventional Fermi- 
liquid corrections as a possible explanation. On the other 
hand, the possibility that this temperature increase is due 
to the presence of a large portion of free localized mag- 
netic moments is probably not very realistic; at metallic 
densities that are 3-^4 times higher than the critical 
density n c , the localized moments should be coupled by 



the RKKY interaction. Here we look at the observed 
temperature dependence from the point of view of the 
non-analytic corrections discussed in the paper. 

The observed trend of the temperature dependence 
corresponds to decrease of the spin susceptibility with 
temperature. Equation ^ predicts, however, an op- 
posite trend. The explanation of the experiment may 
be attributed to the non-monotonic behavior of the spin 
susceptibility due to temperature rcnormalizations of the 
negative amplitudes j n , as has been discussed above. In 
Ref. the data are presented for n/x(T), where n is 
the density of the 2D electron gas. For silicon MOSFET 
(there are two valleys) the density n = p%/n. For the de- 
generate Fermi gas with two valleys the unrenormalized 
spin susceptibility xo = 4f and the ratio n/xo — Pf/^ 171 - 
In the presence of two valleys the expression for 5x as 
given by Eq. JT]) should be multiplied by the factor 4, 
i.e., Sx = 4^(T/ej?)7 C (T) 2 . Since the temperature cor- 
rections are relatively small, one can expand n/x- 

S(n/ X )n^(-5 X ) = -T 1 C (T) 2 . (5) 

The main feature of this relation is that the temperature 
dependence of 6 (n/x) is determined by the amplitude 
7 (T) only. Notice that the modification of the spin sus- 
ceptibility by the Stoner factor drops out from 5 (n/x) 
because of the renormalization of the two external ver- 
tices in Sx, which are indicated by triangles in Fig. [TJ see 
Eq. (lB9ll . 

The data— exhibit a noticeable increase of 6 (n/x) 
starting from T > 2K and correspond to |7 C (T)| ~ 1 
in Eq. (J4|. Such a value of the amplitude 7 C (T) in 
the Cooper channel implies that the experiment has 
been performed at temperatures comparable with the 
temperature of the superconducting instability. Note, 
however, that in the disordered system studied in 
Ref. [l(| the superconducting instability cannot fully de- 
velop for a non-zero harmonic 7n(T) = l/[(l/r^')~ 1 + 
ln(A/ maxT, T~, ast ■ )], because the instability is blocked 
by disorder when r~c t . > T c ; here T c is the temperature 
of the Cooper instability in the dominant harmonic. 

The experimental curves well in the metallic region 
correspond to ep ~ 30 -f- 40 A'. We evaluate the bare 
value of the amplitude j c as r c ~ —0.25-=- -0.3. Then, 
7 C '(T) reaches the value 7 C (T) = 0.5, where Sx changes 
its sign, at temperatures ~ 10AT. Would it not be 
blocked by disorder, the superconducting instability may 
occur at T c ~ 1AT. The measurements are performed 
up to T < AK. In the suggested explanation, the ini- 
tial decrease with temperature of the spin susceptibility 
should be succeeded by an increase at larger tempera- 
tures, T > 10AT. The temperature range of the existing 
measurement is too small, however, to make this conclu- 
sions definite. 

Two comments may be in place here, (i) The analy- 
sis of Sx(T) basing on Eq. @ is only qualitative when 
|7 C (T)| ~ 1. In fact, the situation is more complicated. 
The derivative dx/dT is a series in ji(T) which does not 
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reduce to Eq. ((4]): 

dSx/dT = (u/e F )J2 [{-l) l ' +l "lv{T) 1V .{T) 



ll ,{T) ll „{T) ll ,„{T) + - 



(6) 



The details arc in the end of Appendix [Cl To get a final 
conclusion concerning a non-monotonic behavior, one has 
to calculate for |7 C (r)| ~ 1 the whole series in Eq. (|6|), 
which is hardly possible, (ii) Initial values of the non- 
zero harmonics in the Cooper channel are highly non 
universal. In addition, the destructive influence of the 
disorder can stop the renormalization of the amplitude 
of the dominant harmonic before it becomes noticeable. 
Taken together, these facts mean that the experimental 
situation may strongly vary from sample to sample. 

In a separate publication^, we present an alternative 
explanation of the experimental situation. We consider 
repulsive scattering amplitudes only, > 0, that is, per- 
haps, more appropriate for a system studied in Ref. [10| . 
Since the repulsive scattering amplitudes scale to zero, 
a more delicate analysis is necessary there. One has to 
account for all three channels, i.e. to include the rescat- 
tcring of particle-hole pairs in the zero-sound channel 
in addition to the Cooper and 2pj?-scattering channels. 
In this way it is possible to explain the observed tem- 
perature dependence of the spin susceptibility, both in 
sign and temperature. (Note that the non-monotonic 
behavior may also occur here. The temperature depen- 
dence changes its sign when the logarithmic suppression 
of the repulsive amplitudes become ineffective, i.e., when 
l/ln(e F /T)> 7 c .) 

IV. CONCLUDING REMARKS. 

A decade ago, the authors of Ref. Q made an im- 
portant conjecture, based on the power counting, that 
thermodynamic quantities in the 2D electron gas contain 
linear in T dependences. They argued that similarly to 
one dimension, where the spin susceptibility depends log- 
arithmically on the temperature,— the anomalous tem- 
perature terms should also exist in higher dimensions, 
albeit in a weaker form. 

We believe that presented calculation of the spin sus- 
ceptibility in 2D does not support this line of reason- 
ing. Indeed, the temperature dependence of the spin 
susceptibility in one dimension (ID) originates from the 
logarithmic in temperature dependence of the Stoner 
enhancement calculated in the parquet renormalization 
scheme! 23 ' 24 In contrast to the Fermi liquid description 
(which is analyzed in terms of a propagation of two quasi- 
particles) and the parquet approximation (which is an- 
alyzed with the Sudakov cross-section of two Green's 
functions), the anomalous terms in the thermodynamic 
potential in 2D develops only when two (or more) edge 
parts of the two-particle correlation functions act in the 
combined way. In particular, in 2D, unlike ID, the tem- 
perature dependence of the spin susceptibility does not 



reduce to a mere renormalization of the Stoner factor. In 
our opinion, the physics of the temperature corrections 
in 2D is rather unique and is not a mere continuation of 
ID logarithms. 

It has been shown in Appendices |A1 and [Dl that both a 
2pp-pair and a two-particle section in the Cooper channel 
have non-analytic parts (edge parts). In two dimensions 
the square root singularities of these edge parts are rel- 
atively weak. However, taking derivatives with respect 
to an external parameter, a magnetic field in our case, 
makes them more singular, i.e., activates them. Because 
of high sensitivity to the temperature and magnetic field, 
diagrams containing edge parts running on each other 
produce the anomalous terms in the spin susceptibility. 
When a diagram for the thermodynamic potential con- 
sists only of two pairs of quasiparticles, their edge parts 
are effective because the quasiparticles are constrained 
by the conservation of the momentum and energy. In 
this case the edge parts of the two pairs act together 
in a combined way. For more complicated diagrams the 
non-analytic pieces may be not constrained. Then they 
are smeared out by independent integrations and become 
ineffective. This leads us to the conclusion that for the 
anomalous temperature dependences the essential dia- 
grams in the thermodynamic potential are ladder dia- 
grams only. For the spin susceptibility the set of essential 
diagrams in the 2pp and Cooper channels is organized in 
a way that the anomalous terms can be considered within 
the Cooper channel only. (This rearranging is necessary 
to avoid the double-counting of the two-section terms 
which, as has been found here, completely overlap in the 
2pp and Cooper channels.) We show next that one has 
to keep arbitrary number of sections in the Cooper chan- 
nel ladder to collect all non-analytic terms in the spin 
susceptibility, i.e., not only a product of the two pairs. 
(This is because the singular part of a two-particle sec- 
tion in the Cooper channel depends on the ratio uj / qvF-) 
The truncation of the series in powers of the renormal- 
ized amplitudes is possible only when these amplitudes 
are small. 

In Refs. 013 the possibility of strong renormalizations 
in the Cooper channel was neglected and the ampli- 
tude r(7r) was analyzed from the point of view of the 
Fermi-liquid renormalizations. These renormalizations 
are not related to the infrared logarithms, but they may 
be large near Pomeranchuk's instability in one of the 
channels, e.g., near the Stoner instability in the case of 
the magnetism.— Since the amplitude r(7r) involves all 
angular harmonics, an instability in any of them can be 
important. In this paper we assumed that the discussed 
Fermi liquid system is not in the immediate proximity of 
any instability. In this respect our analysis emphasizing 
only the infrared renormalizations in the Cooper channel 
may be incomplete. 

Finally, let us point out that when the Fourier har- 
monics amplitudes in the Cooper channel, T^, arc re- 
pulsive and scale to zero at low temperatures, the con- 
tributions to the spin susceptibility from the Cooper 
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channel become unimportant. This situation is analyzed 
in a separate publication^, which is complementary to 
the present one. For vanishing Cooper amplitudes, the 
anomalous temperature terms are dominated by the non 
analytic contributions from the particle-hole pairs with 
small momentum transfer (zero-sound channel). Extend- 
ing our analysis to particle-hole rescattcring, we show 
that, exactly like in the Cooper channel, the two-section 
term in the zero-sound channel is dominated by the back- 
ward scattering. To avoid the double counting which 
now involves all three channels, the two-section term has 
been fully attributed to the Cooper channel, where it gets 
killed off by the logarithmic renormalizations of the re- 
pulsive amplitudes. We thus conclude that for the repul- 
sive interaction the anomalous temperature corrections 
are determined by the ladder diagrams in the zero-sound 
channel with three and more rescattering sections.— 



APPENDIX A: ANALYTIC PROPERTIES OF 
THE PRODUCT OF TWO GREEN'S FUNCTION 
WITH q^2p F 

Consider a product of two Green's functions summed 
over the momentum and the energy as a function of the 
transferred energy iuj n and momentum q 



[GG]i 



d 2 p n F (e p+<i - p) - n F (e p - p) 
(2n) 2 e p+q - e p - iw n 



(Al) 

where n F is the Fermi-Dirac distribution function 
n F (c) = [cxp(/3e) + We are interested in the an- 

alytic structure of Eq. (|A1|) in the complcx-o; plane when 
q = |q| f=a 2pp. To get an idea about the analytic proper- 
ties of a 2pi7--pair of Green's functions, let us start from 



the case when e r 



(pj+p±)/2m. We 



first shift the mo- 



mentum p + q to p in the first term in the numerator 
and write the dp-integration as d 2 p = dpudp_i_, where py 
is along the direction of q. Then, the variable p± drops 
from the denominator and enters only via n F (e p — p). 
This paramctrization allows us to introduce the "spec- 
tral density" p(j>\\) which determines the dependence on 
the chemical potential and the temperature 



Pip 



dpi 
9 

-r*r\ " 



n F {e p - p) 



(A2) 



After that get an expression for the 2pp-pair in the 
form of the Lehmann-likc representation in the complex- 
es plane 

[GGUv = — (n,(g >W ) - Ii r {q,oj)) , (A3) 

where the functions Hi{q, u>) and H r (q,to) arc defined as 

m 1 / 2 „ , m [ dp\\ p(p\\) 

p\\ — (m/q)(u! ± q 2 /2m) 

(A4) 
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FIG. 2: Illustration of the analytic structure of the 2pr pair of 
Green's functions. The horizontal axis is Re ui and the vertical 
axis is q. The parabolas show the position of the edges of the 
branch cuts of II;. r in the q, Reoj-plane. The dashed lines show 
the branch cuts of n r .; for two typical values of q: the II r -part 
has branch cuts between the two right-looking parabolas, and 
11; has branch cuts between the two left-looking parabolas. 
The rectangle marks the region where the edges of the two 
branch cuts of two 2p F pairs are close to each other. The 
meeting of two edges leads to non-analytic behavior of the 
thermodynamic potential. 



Here ± in the denominator corresponds to 11/ and n r , 
respectively. From now on, in this Appendix we consider 
the temperature to be zero. [The non-zero temperature 
can be restored with the help of equation (|B3[) as it is 
explained below.] At T = the integral in Eq. (|A2|) is 

easily evaluated: p(p\\)t=o = \Jpp ~ pf^pp— p\\). In the 

complex-w, plane H^ r (tu) considered as a function of the 
variable u± = {raj 'qpp)(tu ± q 2 /2m), have the following 
analytic form 



,1/2 



Attvf 



mp F 

II;. r = F{u±) 



F(u) 



2nq 
1 dx y/1 — X 2 

.1 7T X — U 



—u + \/ u — rv u + 1 . 

(A5) 



Here the square roots should be understood as having a 
branch cut directed to the right on the real axis. Corre- 
spondingly, F(u) has a branch cut between —1 and 1 in 
the complex-it plane. On the upper side of the branch 
cut the imaginary part of F(u) is positive; the real part 
originating in the square root term is positive to the right 
of the branch cut and is negative to the left of it. 

Only the second term in F(u) containing the square 
roots can be responsible for non-analytic behavior in the 
thermodynamic potential. The non-analytic temperature 
dependence in the product of four Green's functions orig- 
inates from region marked by dashcd-line rectangle in 
Fig. [2] where the edges of the branch cuts of n r and n; 
are close to each other. This makes the contribution from 
this region very sensitive to the magnetic field and tem- 
perature smearing. 

For the calculation of the non-analytic term in the spin 
susceptibility it is enough to use an appropriate approxi- 
mation for 2pi?-pair in the region marked by a dashcd-line 
rectangle in the Fig. [5] This involves the u = — 1 edge of 
branch cut in n r and the u = 1 edge of the branch cut in 
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IX; . In the following we omit the smooth terms in F(u). 
It can be checked by an explicit calculation that the first 
term in F(u) docs not to contribute to the anomalous 
temperature behavior. 

The spectrum of electrons is not considered to be 
quadratic, like in Eqs. (|A4j) and (|A5|) . but we still as- 
sume the spherical symmetry. The generalization for a 
non-spherical spectrum can be done, if needed. To devise 
an adequate approximation for II/ r , we consider again 
the product of two Green's functions, [G + G_] g de- 
scribing the propagation of two quasiparticles with the 
momenta p + and p , which are almost opposite to each 
other. The transferred momentum q = p + — p_ is close 
to 2pp, and 5q — |q| — 2pp describes a small devia- 
tion from 2pp. We choose the following parametrization 
where pj^ is a vector directed along q: 



q 
p 



2p F - 
Pf 



<5q, 
- Sp . 



(A6) 



To reproduce the correct threshold behavior we keep min- 
imal terms in the expansion in pii and p±. We linearize 
the energy spectrum around the points ±pp in the direc- 
tion perpendicular to the Fermi surface, but preserve the 
quadratic momentum dependence in the direction paral- 
lel to the Fermi surface to account for the Fermi surface 
curvature 



describes the square root function with the branch cut 
directed to the left on the real axes of w-plane, while n r 
describes the function with the branch cut directed to the 
right. [Indices / and r correspond also to the direction 
of "propagation" in the combinations uj ± Vpq.] Both II; 
and II r have a positive imaginary part on the upper side 
of their branch cuts in the complex cj-planc. 

Let us conclude this Appendix by a comment on 
Eq. (|A9j) . This equation gives a representation of the 
edge parts of 2j>F-pair of Green's functions as two pieces, 
n r and III , in the region inside the dashed-line rectangle 
in Fig. [2] Together these edge parts act as a sort of an 
"edge-mode" ; we mean the combination uj ± vp(q — 2pp) 
in the denominators of Eq. (|A9[) . The square root spec- 
tral weight of these edge-modes is much weaker compared 
to the usual quasiparticles, e.g., phonons. However, tak- 
ing derivative with respect to an external parameter, a 
magnetic field in our case, makes the spectral weight sin- 
gular. As a result, these edge modes become active. In 
Appendices IBI and ICl we show that, because of the un- 
usual sensitivity to the temperature and magnetic field, a 
diagram in the thermodynamic potential with two propa- 
gators of the edge modes running on each other produces 
the anomalous temperature dependence in the spin sus- 
ceptibility of the interacting electron gas. 



■ p+ 



fi = - v F Sp\\ +p 2 L /2m 

Ho = v F (5p\\ + 6q) + p\/2m. 



(A7) 



We obtain the following approximation for H r i (ui) at zero 
temperature 



,1/2 



4irvp 



n;, r (g, u) 



1 

2^ 



dp\\dp ± np(±v F (Spi\ +6q)+p\/2m-5fi) 



2- 



vf(25p\\ + Sq) — lo 



(A8) 



where ± corresponds to IT; and to II r , respectively. 

Having in mind the calculation of the spin susceptibil- 
ity from the correction to the thermodynamic potential 
we introduce a shift of the chemical potential fi = no + d/i 
in Eq. (|A8j) . Introducing x = 2vf(Sp\\ + Sq) ± 26 /i in H r 
and III , respectively, we obtain after integrating over p± 
the following expressions for n/. r 



Ui(6q, lo)t=o 



dx 



U r (Sq,uj)T=o = J ~ 



7r x — (o> + VF&q — 25fi) 
dx ^fxQ(x) 



x — [ui — VF&q + 25fi) 



(A9) 



These expressions correspond to Eq. (|A4[) within the re- 
gion marked by dashed-line rectangle in Fig. [5J Here the 
functions H;. r arc defined with positive square roots in 
the numerator. In this Lchmann-like representation, n; 



APPENDIX B: ANOMALOUS TERMS 
GENERATED BY TWO 2p F -PAIRS. 



The term in the thermodynamic potential with four 
Green's functions which produces the anomalous term in 
the spin susceptibility is 



Jfi 2 (T, h) = -(1/2) S ^(T; h, -h)\I(n)\ 2 . 



(Bl) 



where g^Sfjti, 5/j,2, T) signifies the product of the two 
pairs of Green's functions with the chemical potential 
shifted in two pairs by 6fii and 5fi2, correspondingly; 
arrows indicate spin and we consider a Zeeman energy 
h = (gfiB/2)H as the source of a shift of the chemi- 
cal potentials. The factor 1/2 can be obtained from the 
comparison with the second order term in the perturba- 
tion theory for the thermodynamic potential. [There are 
two such contributions, with fj, and |f arrangements of 
spins in the two pairs of Green's functions, which has 
been accounted for by the factor 1/2 in Eq. (|Blj) ]. The 
terms with the parallel spin projections in the two pairs 
of Green's functions, i.e., g± (T; h, h) and g^(T; h, h), do 
not contribute to the anomalous term. 

Having in mind the correction to the thermodynamic 
potential 5fl2 given by Eq. (jBlj) , we now perform the sum 
over Matsubara frequency iuj n . The standard manipula- 
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tion leads to 



g±{T; £/zi,<5// 2 ) 



mpp 



dSq f duj (3oj 
7T / 2?r 2 



(4ttv f ) 2 

{im [n,(^i) - n r (^i)] T Re [n,(fy a ) - n r (^ 2 )] T 

+ (fyi ty 2 )} • (B2) 

The signs of different terms in this expression are com- 
pletely determined by the analytical properties of func- 
tions H rj i- Only the cross-product terms, i.e., II; n r , be- 
have sharply in the vicinity of u) = 0, and therefore are 
responsible for non-analytic terms in the thermodynamic 
potential. It can be checked directly that in a product 
of, let us say, II r with a regular function its square root 
singularity disappears after the ^-integration. For sim- 
ilar reasons, the products II; II; and IL,II r yield regular 
functions in frequency. On the contrary, in the product 
II; n r the edges of the two branch cuts run on each other, 
preventing the smearing of the singularity. 

The temperature dependence in the functions II; r can 
be restored with the help of the relation 



n(w) T 



ds - 



driF(s) N 



n(w - s) T=0 (B3) 



which follows from the identity 



n F (co) T = J ds rip {lo - S ) T=0 . (B4) 



Performing the g-integration in Eq. (|B2() one gets 

fx x \ 1/3 (a a a dn F ( Sl ) dn F {s 2 ) 

g4{dfii,dfi 2 )T = -o — / dudsids 2 

o€f J 



coth — 
2 



there are no "drag" diagrams for the transverse spin sus- 
ceptibility, because Green's functions attached to the ex- 
ternal vertices a ± , carry opposite spins. The spin suscep- 
tibility correction which corresponds to these diagrams is 



Sx = 2( X oM 



ld 2 g 



d8n\ 



Sfi 2 =Q 



TO 



(B6) 

Here the ratio {xo/ v ) accounts for the factors arising in 
the calculation of the spin susceptibility from the corre- 
sponding diagram with cr + and a~ at the vertices. Factor 
2 accounts for a spin trace in the upper bubble in the dia- 
grams presented in Fig. [1] The derivative in the brackets 
generates an expression with six Green's functions corre- 
sponding with the correct coefficient. 

The to integration in Eq. (|B6|) can be cast in the form 



1^4(1;^,^; 
2 dSf4 



V 

8e7 



du> 



2 — (3oj coth 



2 



cosh flu — 1 



(B7) 



The first term here is proportional to the ultraviolet cut- 
off. Since it does not depend on the temperature, it will 
be not considered here. [In fact, the traces of this term 
reveal themselves in the T 3 terms which cannot be re- 
duced to the renormalization of the amplitude r(?r) by 
the rescattering in the Cooper channel.] Omitting the 
factor (g/is/2) 2 , we obtain for the anomalous contribu- 
tion in the spin susceptibility 



ds i ds2 

(uj + 8^1 — 5^2 + si — s 2 ) 2 sign(cj + A + si — 

(B5) 



Sx = v |^r(7r) 



T 



(B8) 



The obtained expression depends on the relative shift 
A = 5fj,i — 5fj,2 of the chemical potential in the two 
pairs of Green's functions. The signum function in 
this expression appears because the contributions with 
ImII r (<5/ii) RcII;((5^t 2 ) and ReII r ((5/ii) ImII;(5/x 2 ) have 
opposite signs. The expression (|B5|) contains the needed 
non-analytic frequency dependence through the factor 
signfw + <5/ii — <fyt 2 ). This factor compensates the odd- 
in-frequency behavior of coth allowing for the cubic A 2 T 
term in the expansion of (?4(A,T) which produces the 
linear in T term in the spin susceptibility. Notice that 
the expression obtained in Eq. (|B5[) explains how the de- 
pendence on the order of taking the limits H — > 0, T — > 
appears in the thermodynamic potential. At T <C H 
the temperature dependent terms in the thermodynamic 
potential vanishes as T 3 . To calculate the spin suscepti- 
bility one has to take the opposite limit, H <C T. 

Let us first consider a contribution to the transverse 
spin susceptibility, as illustrated in Fig. [TJ Notice that 



S2 ) which coincides with the result of previous calculations, 
see Refs. @H. 

Let us clarify the role of the amplitude r(0) with the 
scattering angle 6 — 6' « in the linear in T term in 
the spin susceptibility. In the diagrams with two pairs 
of Green's functions, which can also be looked at as two 
pairs with the momentum transfer close to 2pp, the scat- 
tering amplitude T(0) can be mixed with r(7r) or used 
twice alone. With one exception, these are one-loop dia- 
grams which cannot contribute to the spin susceptibility 
(all spin projections are the same and they don't have 
part with a difference in the chemical potential). The 
only exception is the diagram with two loops. In each 
loop the two Green's functions have close momenta, but 
their directions arc opposite in two loops; the spin pro- 
jections are also opposite in the two loops. When rear- 
ranged as two pairs with the momentum transfer close 
to 2pp, they have needed shifts of the chemical poten- 
tials, but still 5jii = S/j.2, and therefore the corresponding 
gi(6fJ,i,6fj,2) does not contribute to the spin susceptibil- 
ity. With this comment we fixed the arrangement of the 
dashed lines in Fig. [TJ We see that the magnetic field 
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dependence drops out from unless there are two loops, 
each with a 2p^-pair of Green's function, and with the 
spin projections directed oppositely in each of the loops. 

The performed calculation of the anomalous term in 
the spin susceptibility demonstrates that it is generated 
by the pole parts of Green's functions close the Fermi 
surface. It follows from this fact that the Fermi liquid 
renormalizations of the 8\ originating from the two ex- 
ternal vertices (triangles in the Fig. [I} are equal to the 
square of the Fermi liquid parameter 1/(1 + Go) respon- 
sible for the renormalization of the Pauli spin suscepti- 
bility. Including these renormalizations the temperature 
correction takes the form 



(1 + Go) 



(B9) 



The longitudinal spin susceptibility can also be ob- 
tained in the same way by taking the second derivative 
with respect to the magnetic field, \ = —d 2 Vl2/dh 2 . No- 
tice that in the case of the transverse spin susceptibility 
the coefficient 2(xo/v) originates from two independent 
spin traces in Fig. [TJ In the case of the longitudinal spin 
susceptibility there is no spin traces anymore, but the 
additional factor 4 in this case originates from the dif- 
ferentiation of gi(T\ h, —h) in Eq. (|B1[) with respect to h 
rather than with respect to its first argument only as in 
Eq. (|B6[) . As a result one obtains the same expression as 
in Eq. (jB8|) . 



APPENDIX C: NON-ANALYTIC TERMS IN THE 
COOPER CHANNEL 

In Appendix [B] it has been assumed from the very 
beginning that the transferred momentum q is close to 
2pp. We start the present Appendix with an alternative 
treatment of the term (50 2 which now will be recalcu- 
lated as two sections in the Cooper channel. Here we 
will not assume a priory that only backward scattering 
is important. Remarkably, we obtain this fact as a re- 
sult of independent calculation. Finally, we extend the 
consideration of the anomalous terms by analysis of the 
Cooper ladder in the thermodynamic potential. 



1. Calculation of two pairs in the Cooper channel. 

We start with a pair of Green's functions in the Cooper 
channel 



[G(tyi)G(*A*2)],,to 
d 2 P 



where £ p = e p — fx and q = |q|. To account for the 
effect of magnetic field, we introduced the shift of the 
chemical potential 8^ = ±h depending on the sign of spin 
projection. We now shift the energy variable £ p by (dfi\ + 
<5/Lt2)/2. The obtained expression depends on the relative 
shift A = 8m — S/j,2 of the chemical potential in the two 
Green's functions. Only the pairs with the opposite spin 
projections retain the dependence on the magnetic field 
h; the anomalous temperature terms are not sensitive to 
the overall shift of the chemical potential. Finally, after 
making analytical continuation to the complex-w plane 
iu>„ — > lj wc obtain: 



[G(^i)G(<W; 



d 2 P 



2 J (2n) 2 
A<-> -A. 



2 tanhf (£ + (qv F /2)cosd - A/2 



2£ 



(C2) 



where 8 is the angle between the momenta p and q in 
Green's functions in Eq. (|C1[) ; here we expand £ p + q /2 
using the smallncss of q and introduce £ = £ p ; i.e., 
Cp+q/2 = £ + (qvF /2) cosO. The result reproduces the 
Cooper logarithm multiplied by the single particle den- 
sity of states (per spin). We use a variable x — 2£ to 
obtain the expression in form of the Lehmann-like repre- 
sentation: 



[GG] q 



where 



^ 9,u + A) T + *(g, 9,u- A) T 



(C3) 



*(g,0,w± A) T = 



dx tanh j (x ± A — qvp cos 0) 



_(C4) 

The normalization of ^(0) is such that the imaginary 
part of \P approaches ±1 at large values of u. At zero 
temperature the hyperbolic tangent becomes the signum 
function. The temperature width of Fcrmi-Dirac function 
can be restored from T = function with the help of the 
relation 



8 tanh £2 , 
ds — (oj - s) T =o 

OS 



that follows from the identity 
tanh ax 



1 f d tanh as 
2/ dS ^^ Slgn(x - s) 



(C5) 



(C6) 



The contribution to the thermodynamic potential from 
two sections (recall that by section we understand a pair 
of Green's function describing a propagation of two quasi- 
particles between the scattering events) in the Cooper 
channel is given by 



1 f d 2 p tanh | (£ P + q /2 - 8m) + tanh § (£- p+q /2 - 8fi 2 ) 

2 J (2?r) 2 ?p+q/2 + £-p+q/2 - ^Ml ~ 8fJ, 2 - i^n 

(CI) 



n 2 (A) = -i(f 



qdq d#i c^2 j,y 
2tt 2tt 2tt f-* 



e 2 )r( 



) 2 ,q 1 iU} n )T 



Oi)T( 



'l,q, iVn)T , 

(C7) 
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where T(A) T = + A) T + - A) T ]/2, and 

r (0) = ^r(0) is dimcnsionless angle-dependent ampli- 
tude in the Cooper channel. Let us pass to the Fourier 
harmonics of the Cooper channel amplitudes, T c (0) = 
J2 n e mS T^. We assume that the amplitudes T(9) are real 
and are even functions of 81 — 62- This corresponds to real 
r' = r~'. The angular integration leads to the following 
expression for the contribution of two Cooper sections to 
the anomalous part of the thermodynamic potential 



1 



qdq T V 



T; Ti_ m (g,io;„)Tr m T m _i((j, iu n ) 



t ■ 



(C8) 



Here the angular harmonic ^> n of function \&(#) are de- 
fined as 



(C9) 



After transforming the frequency sum to the integral we 
obtain 



Q 2 (A) = -£^ 



dio coth —— 



Im {T„_ m (<7, u>, A) T T m _„(g, u>, A) T } 



(CIO) 



where the contour of (^-integration is slightly shifted 
above the real axes. 

Let us now turn to the calculation of the anomalous 
temperature term in the spin susceptibility. For the 
chemical potential shifts induced by the magnetic field 
A = 2h and, therefore, 5\ = - Jj^Q, 2 {2h) . This yields 

S X "4^2" J duicoth ^f r n r ?n (h nm (u;s 1 ,s 2 )) T 

+ (W1«W2). (Cll) 

From now on, q is an energy variable, i.e., qvp — > q. This 
substitution leads to the factor l/v F = nv/ep in the first 
line and makes h nm to be dimcnsionless. The notation 
(■■ -}t means the integration over si,2 which originates 
from Eq. (|C5j) . In h nm the shifted frequencies u>± >2 = 
w — si,2 have been introduced and here (and everywhere 
below) *(cj) means *(u>) T =o- Next, in Eq. (|C11|) we 
replace the second derivative with respect to A by the 
second derivative with respect to u> and afterwards the 
limit h — -> has been taken. Notice that terms with the 
first derivatives vanish because dT(h — > 0)/dh = 0. The 
q-intcgrations in h nm yields (for details see Appendix [Dj 



-2sigii(wi + W2)(-l)' 



(C12) 



It is now necessary to perform the u> and s% t2 integrations 



coth 



T%Tg(-l) n+m fd Sl ds 2 



2nv 2 F 



'-duj 



BudtaahPp dtanh^ 

sign(2w - si - s 2 ) . 



dsi 



ds 2 



(C13) 



The w-integration can easily be done with the use of 
the relation coth(/3w/2) = 2Tdlnsinh(/3w/2)/dw. Trans- 
ferring the frequency derivative from lnsinh(/3a>/2) to 
sign(2o! — si — s 2 ) we get a 5-function of frequency. The 
boundary term here corresponds to the first integral in 
the right hand side of Eq. (|B7[) . It does not depend on fre- 
quency and will be dropped. Passing to u\ >2 = tanh/Jsi^ 
as our new integration variables, we get an elementary 
integral. Finally, 



5x = ,-ErK(-i 

€p — 



(C14) 



The angular integration Eq. (|C7[) covers all momentum 
directions in both Cooper sections. Nevertheless, the re- 
sult given in Eq. (|CJ14|) can be written in terms of the 
backward scattering amplitude only 



6 x = v-\r c (n) 



(C15) 



where T c (ir) = V (-l) n r^. It is crucial for this result 
that, remarkably, the q- integration in h nm depends on 
harmonic indices only through the factor (— l) n_m . This 
also implies that the scattering amplitude T(0) drops out 
from the anomalous temperature corrections to the spin 
susceptibility calculated to the second order in Tf, see 
Ref. [y. The result obtained here coincides with Eq. (|B8|) . 
which has been calculated within the 2pp channel. 



2. Cooper ladder in the thermodynamic potential. 

The anomalous terms in the spin susceptibility gener- 
ated in the Cooper channel are described by the following 
expression 

Sx = J^ d0J coth Tr Im i d ^ {q ' Ld)T ^ c ) ' 

(cie) 

where 1 



7c = L c [l + (7r/2)§(g,u>) T fc] 



(C17) 



represents a matrix of harmonics 



r^fSrim is a matrix of the Cooper 



Here (*)„ m = 
of $(0) and (f c ), 
channel amplitudes. The trace is over angular harmonic 
index n. Notice that only the terms with the second 
derivative with respect to magnetic field acting on the 
same Cooper section survive the limit h — > 0. 
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Let us first explain why the term d^o(uj) alone does 
not contribute to 5\ (Obviously, such term cannot ex- 
ist for non-zero harmonics). Using relation {u>d q + 
<ZcL)<9 w <Fo(<7, w ) = one can reduce the g-integration 
to the boundary terms which however vanish because of 
absence of the imaginary part. Therefore, at least one 
more \1/ is necessary to be a partner of d^,^o(q, w )- After 
choosing the partner in the matrix 7c the expression in 
Eq. ([TO]) takes the form: 

x frtj , u fa 
ox = / — — old coth — — 
/ 2tt 2 



the integrals which determine this quantity converge at 
lu ~ T and depend on the ultraviolet cutoff only through 
the renormalizations (compare with Eq. (|B7j) ). The dif- 
ferentiation of the Eq. (|C19[) gives 

dSx/dT = (^/e F )^(-l) n+m [7n(T) 7m (T) + 2 7 „(T) 2 7m (T)] 

nm 

(C23) 

In fact, the situation is more complicated. The deriva- 
tive dx/dT is a series in 7;(T) which does not reduce to 
Eq. (IU231) : 



x Trim 



{dl*{q,w) T %*{q,w) T %) . (C18) dSx/dT = (l//ejp) £ [(-!)''+<%, (T) 7 *»(T) 



Now let us point out the special role of the zero har- 
monic in 7,7. Any {^jc)nm can be constructed starting 
from r^\I> n _;< • • • ^///_ m r^ which is chosen to be irre- 
ducible with respect to $f . Each F^ can be dressed 
by arbitrary number of "Jo by replacing with 7„ = 
r£/(l + (7r/2)r^* ). The specifics of * is that it de- 
pends logarithmically on the ultraviolet cutoff. Let us 
single out this term, i.e., "Jo = ln(A/T) + ^ , o- As we know 
from the experience of the second order calculation the 
anomalous contribution accumulates from uj,qvF ~ T. 
Therefore <Fo ~ 1 a s well as "F„^ m . Now neglecting 
'Fo and "Fn^m in 7c we come to the expression for the 
anomalous term in the spin-susceptibility, which is renor- 
malized by Cooper logarithms 



T 



^(-l)"7r,(T) 



where the renormalized "f n (T) is 



7n(T) 



1 + r^lnA/T' 



(C19) 



(C20) 



This result is in full correspondence with the renormaliza- 
tion group equation for the angle-dependent amplitude 
T c (9) in the Cooper-channel 



dT c {6) 



dO 1 
2^ 



T C {0-6')T C {6') 



(C21) 



where £ = InA/T. [Notice that apart from renormaliza- 
tions of r^, the large logarithmic part cannot contribute 
when a partner of in Eq. (|C18[) is chosen to be the 
zero harmonic, because of the same argument as pre- 
sented in the beginning of this paragraph.] 

The result given in Eq. (|C19[) can be presented in terms 
of the renormalized backward scattering amplitude 



T 



5 x = v—\r c (n) T \ 2 



(C22) 



where T c (tt) t = Yj n i^ l Tln{T). For a general analy- 
sis of the temperature dependence of 5x(T) it is more 
appropriate to work with dSx/dT rather than with the 
spin susceptibility itself. Unlike the spin susceptibility, 



lv {T) llll {T) ll „,{T) + 



(C24) 



This series is generated when in Eq. (|C16[) an arbitrary 
number of ^ n ^m or 'Fo are taken as a partner for <9 2, F, 
and it differs from Eq. (|C23|) . The fact that they are 
different is obvious when all VI" ■ ■ ■ are not equal to each 
other, but this is also valid when some (or all) of Vs are 
equal. The reason why the series in Eq. (|C24[) is more 
sophisticated than a mere differentiation of Eq. (|C19[) is 
because the functions ^ n ^m and 'Fo depend strongly on 
the ratio uj j qvF- In the microscopic Fermi liquid theory^ 
the scattering amplitude in the zero sound channel has 
exactly the same feature. 

Let us point out to a rather subtle contribution to C3 
which gives additional evidence that a simple renorm- 
group generalization of the second order result given by 
Eq. (|C19[) is not complete. As it was indicated in Ap- 
pendix [B] the calculation of the product of four Green's 
functions contains the term with the integration which is 
limited by the ultraviolet cutoff. This term did not con- 
tain any temperature dependence and has been dropped. 
However, when matrix elements of 7c (u>) acquire fre- 
quency dispersion (see Eq. (|C17[) ) the result of this inte- 
gration ceases to be a featureless constant and it produces 
a new contribution to C3, which does not reduce to the 
renormalization group generalization of the second order 
result Eq. ([CT9|) . 



APPENDIX D: CALCULATION OF THE 
FUNCTION h nm (ui,wa). 

The function 'Fo(w) is rather special because it con- 
tains the logarithm of the ultraviolet cutoff. Therefore, 
we will consider the contributions involving zero and non- 
zero harmonics separately. We show that the result of the 
g-integration in h nm depends on harmonic indices only 
through the factor (— l) n ~ m . This remarkable feature 
is responsible for the fact that the linear in tempera- 
ture term in the spin susceptibility calculated from two 
Cooper sections depends only on the backward scattering 
amplitude r c (7r) = £(-l) n r£. 
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1. Terms with zero harmonic of 9(9). 

To obtain an expression for zero angular harmonic 
we first calculate its imaginary part at frequencies shifted 
slightly above the real axis. At zero temperature we get 
from Eq. ([C4]l for lu = u + iS: 



M > q 

\oj\<q (Dl) 



with the w-derivativc. The result is 



Im ^(g, oj) = signer, 

= (2/ir) arcsin(w/<?), 



The corresponding analytic expression for ^>o applicable 
in the complex-cj plane is 



*o(9,w) 
= { — in + 

7T I 



In 



+ (m 2 - q 2 ) 1 ' 2 
A 

( ^uu + {uu 2 -q 2 ) 1 ' 2 



A 



} (D2) 



where A is an ultraviolet cut-off. The notation ln^^z) 
and ln'*~'(z) signifies a particular choice of the com- 
plex logarithmic function. On the real axis the function 
ln^-'(z) has a branch cut at z < and it is real at z > 0. 
The function ln^ - ^ has a branch cut on the real axis at 
z > 0; its imaginary part is equal to zero on the upper 
side of the branch cut and to 2iri on the lower side of 
the branch cut. Altogether, near the real axis of w the 
function Im "to = isigno; for lu = lu ± S when \tu\ > q. 
The square root y 'us 2 — q 2 in the complcx-tj plane has a 
branch cut between q and -q on the real axes; it is real 
and positive at to > q and negative at lu < —q. 

Let us consider the diagonal part h nn . 

/i„„(w) = | J qdqlm{ d^ (q,ui) *o(g,w 2 )} 

+ (lUl <-> LU 2 ) 

= ~\ J dquJl § q /[m { d ^o(q,^i) *o(<7,^2)| 

+ ^ J dglmj d u! i b (q,uj 1 ) uj 1 -^^ (q, w 2 )| 

+ (u>i «-> lu 2 ) ■ (D3) 



To eliminate the second derivative of (^i ) we used the 
relation (qd u + Lud q )[d^ a ] = for [9 W * ] = c/{lu 2 - 
q 2 ) 1 ! 2 to replace the lu derivative with the g-derivative. 
From now on, we use c = —2/tt. In the right-hand side of 
Eq. (|D3|) the first term is constructed to be a full deriva- 
tive. 

Now, using the relation [uid^ + qd q )^o — c which fol- 
lows from the explicit form of ^! , we replace in the last 
term on the right-hand side of Eq. (|D3|) the q-dcrivativc 



h nn {uj) = -{ 



uJidq-^-lm ^ duj^Q(q,uJi) # (<7,w 2 )} 



+ |y dqIm{[d u *o(Q,"i)]^y)} 

+ (u>l <-> LU 2 ) ■ 



0J\UJ 2 

q 



d^o(q,ui 2 ) } 
(D4) 



The integral of the full derivative in the first term is non- 
zero at the lower limit. Having in mind ^[d^^q — > 
0, lu)] — c and Im ty(q — > 0, lu) = signw we obtain for it 



2 ' 

signw 2 ■ 



dq §q lm { dw ^°^ Ul ) *o((?,W2)} 



(D5) 



To calculate the integral in the second line we consider a 
function c/ (lu 2 — q 2 ) 1 ^ 2 in the complex-g plane. For posi- 
tive lu we define the square root in the complcx-g plane as 
a function which has a branch-cut outside [—lu, lu] and has 
a negative imaginary part above a brunch-cut at q > lu. 
With this, we can relate this expression to the integral in 
the complex-g plane along the contour C which envelopes 
the branch-cut in the counter-clockwise fashion. The in- 
tegral is calculated by deforming the contour so that it 
encloses the pole at q = 0. For positive lu\ it is equal to 
— 1; for lui < it changes the sign. Overall, we have 



J dglm{[0 w *o(g,wi)](^)} 



7T 

2 



= ~2 [UJl 



d u ^ (q,uji] 



-signal 



(D6) 



The integral in the last line of Eq. (|D3|) can be cal- 
culated by similar method. A delicate point here is to 
notice all the changes of signs 



LU\LU 2 \ 



q 



J [d^ (q^2)]} 
(D7) 



dqlm | ^* (g,^i) 

= signal + w 2 )sign(u;i)sign(a; 2 ) 

Collecting all terms, we find for h nn 

h n „ =2[sign(wi + cj 2 )sign(cji)sign(w 2 ) - signwi - signup] 
= - 2sign(wi + uu 2 ) (D8) 

2. Terms involving non-zero harmonics of 9(6). 

To obtain a zero-temperature expression for the 
non-zero angular harmonics of ^f(6) we first evaluate 
the angular integral lm iff n (q, to) = J(d6'/27r)sign(a; — 
qvF cos9)e me 

Im*„(g,w) =0, H>q 

= — (2/ nir) sin (n arccos(w / q) ~\ , \u>\ < q . 

(D9) 
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In the complex-^ plane the analytic expression for ^> r 
has a form 



2 (u-(uj 2 -q 2 ) 



w\n\ 



2 _ „2U/2 \ I™ 



q 



(D10) 



We need to calculate the q- integral in h nm (w) 
ftn-m(wi,w 2 ) + /t n _ TO (w 2 ,cJx) where 



/i;(wi,cj 2 ) = - 



gtiqlm | 



if 

9w 



)}• 

0)11) 

Let us sketch the way we perform the g-integration in 
this case. When both u)\ and w 2 are positive, the q- 
integral can be done using contours in the complex-g 
plane. It yields hi(u>i > 0,lu 2 > 0) = -(-1)'. Let 
us look on this result in more details. The integral can 
be written as (tt/2) J qdq[lm ARcB + Re Aim B], where 
A = d 2 ^i{q,uji) and B = *;(g,w 2 )- When w\ < u 2 the 
lower limit for the term with Rc^lImS is q = u; 2 . The 
integral (tt/2) J w qdq(Rc Aim B) is purely trigonometric. 
It is identified as —1 for even I and it is zero otherwise. 
Correspondingly, the integral (tt/2) qdq(lmAReB) is 
1 for odd I and zero otherwise. Alternatively, for w 2 < 
the integral (tt/2) qdq(lm ARe B) is purely trigono- 
metric. It is identified as —1 for even I and zero other- 
wise. We now use this information to resolve the delicate 
question of the sign changes in various domains in u>x, oj 2 
plane. 

A peculiar feature of the g-integration is that at any 
value of ui we take the value of ^ on the upper side 
of its branch cut in the complex-w plane. From the 
definition of \J uj 2 — q 2 given after Eq. (ID2I1 



it fol- 
lows that for u> slightly above the real axis ^>i(—uj) = 
(—l) l ^l(w); the same holds for Ai(lu) and Bi(lu). We 



start with the case when wi < and \lui\ < w 2 . In 
this case the (/-integral in hi(u>i,u>2) has to be defined 
as (tt/2) J qdq(-l) l [-lmARcB + RcAlmB}. Using the 
values for the two g-integrals found above, we obtain for 
this combination again —(—1)'. Therefore, the q-integral 
in hi(u)\,u>2) does not change when w\ changes its sign. 
When |u; 2 | < oj± and u>2 < the q-integral in /i;(a>i,a; 2 ) 
becomes (tt/2) J qdq(-l) l [lmARe B - Re AlmB}. Here 
again, this combination reproduces —(—1)'. Hence, the 
integral remains constant when it is extended across the 
boundaries of the quadrant Wi,W2 > 0. On the other 
hand, when both ui\ and w 2 are negative the integral has 
opposite sign, hi(uJi < 0,u>2 < 0) = (— 1) ; . Altogether, 
it follows from this analysis that hi(uj\, 0J2) changes sign 
when u>i + w 2 = 0; i.e., 

/ 1 „ m = -2sign( Wl +^ 2 )(-l)"- m . (D12) 
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